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ABSTRACT 



Context. The hot subdwarf B star HS 0702+6043 is known as a large-amplitude, short-period p-mode pulsator of the EC 14026 type. Its 
atmospheric parameters place it at the common boundary between the empirical instability regions of the EC 14026 variables and the typically 
j^-^ ' cooler long-period g-mode pulsators of the PG 1716 kind. 

, Aims. We analyse and interpret the photometric variability of HS 0702+6043 in order to explore its asteroseismological potential. 
Methods. We report on rapid wide band CCD photometric observations to follow up on and confirm the serendipitous discovery of 
multiperiodic long-period luminosity variations with typical time scales of ~1 h in HS 0702+6043, in addition to the two previously known 
pulsations at 363 s and 383 s. In particular, we isolate a relatively low-amplitude (~4 mmag), long-period (3538+130 s) light variation. 
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O ' Results. We argue that the most likely origin for this luminosity variation is the presence of an excited g-mode pulsation. If confirmed, 



HS 0702+6043 would constitute a rare addition to the very select class of pulsating stars showing simultaneously parts of their pressure and 



(Z) ■ gravity mode pulsation spectra. The asteroseismological potential of such stars is immense, and HS 0702+6043 thus becomes a target of choice 



for future investigations. While our discovery appears consistent with the location of HS 0702+6043 at the common boundary between the two 
families of pulsating sdB stars, it does challenge theory's current description of stability and driving mechanisms in pulsating B subdwarfs. 
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1. Introduction 

Subdwarf B stars populate the extreme horizontal branch 
(EHB) in the effective temperature range of 22 000 to 40 000 K 
and have surface gravity values from logg=5.0 to 6.2 in cgs 
units. The masses of these hot, evolved objects sho uld cluster 
aroun d 0.5 Mq as suggested by evolution theory jHan et all 
I2003I) . Indeed, masses between 0.457 and 0.49 Mq have been 
derived recently for fou r EC 14026 stars by asteroseismology 
JCharpinet et all 120 05). They are believed to be core helium- 
burning but with hydrogen envelopes too thin to sustain H- 
shell burning. Standard tracks of stellar evolution do not cross 
the EHB region since they do not produce inert hydrogen en- 
velopes, but the high fraction of binaries among the sdB stars 
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* Based on observations collected at the Centro Astronomico 
Hispano Aleman (CAHA) at Calar Alto, operated jointly by the Max- 
Planck Institut fur Astronomie and the Instituto de Astroflsica de 
Andalucfa (CSIC), and at the Bok Telescope on Kitt Peak, operated 
by Steward Observatory (University of Arizona). 



suggests that close binary evolution may play an important role 
in their formation. The relative importance of several proposed 
feeder channels remains unclear. 

To learn more about sdB structure and hence their evolu- 
tionary history, asteroseismology is one of the important me- 
thodical approaches, made possible due to short-period oscil- 
lations exhibited by a fraction of sdB stars (called sdBVs). To 
the group of pulsating sdB stars known as V63 1 Hya variables 
(more commonly referred to as EC 14026 variables) which 
show periods of the order of a few minutes. lGreen et all d2003l) 
have recently added a longer-period group (tentatively called 
IpsdBV) for which the prototype is PG 1716+426. In the short- 
period sdBVs, a k mechanism drives the low-order p-modes, 
where the required opacity bump is due to iron a ccumulated by 
diffusi on dCharpinet et all fl996h . According to lFontaine et all 
(2003J, the same k mechanism can also drive high-order g- 
modes in the long-period sdBVs that have periods of the order 
of about an hour. 

The mfl=15 star HS 0702+6043 was identifie d as a me mber 
of the EC 14026 group in a search program bv iDreizler et all 
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Fig. 1. Original data points (dots joined with a line), and resid- 
uals (thick line) after subtraction of two short-period frequen- 
cies f\ and fi: The variability pattern of the residuals strongly 
resembles the multiperiod brightness changes observed in lps- 
dBV stars. 

Table 1. Frequencies extracted from the February 2004 data. 





frequency [/uHz] 


period [s] 


amplitude [mmi] 


fi 


2754 (9) 


363.1 


21.7(5) 


h 


2606 (9) 


383.7 


4.6(5) 


h 


283 (9) 


3538. 


3.7(5) 



riods (see Tabled were prewhitened from the data. Note that 
the value DR02 have identified for their tentative second short 
period is separated by a one-day alias from the more reliable 
result for fi from the new data. 

Both the full light curve and the residual after subtraction 
of fi and fi are shown in Fig.^ The frequency spectrum of the 
full data is displayed in Fig.|2] A false alarm probability analy- 
sis of the residual (see DR02 for details of our implementation 
of this method) for the peak near 271 juHz in Fig.|2 indicates 
that it is real with a significance well above 5cr. Moreover, 
there is significant residual power left around this peak after 
prewhitening of fy. Regardless of ambiguities in the prewhiten- 
ing due to one- and two-day aliases, this is already evident in 
the original frequency spectrum as a "shoulder" or "hump" to 
the right of the main peak. Since this sub-structure is not re- 
solved from what we have identified as fi, we do not however 
attempt to assign any further frequency detections to it. 

3. What is causing the variation? 

There are three possible interpretations of the low-frequency 
variation: binarity (reflection effect or ellipsoidal variation), a 
rotating star spot, or pulsation. Its coherence and persistence 
over four years rules out normal sky or transparency variations. 



(2002, hereafter DR02). Its effective temperature was deter- 
mined to be 28 400 K, its surface gravity log §=5.35. These 
spectroscopic parameters place it at the cool end of the 
EC 14026 instability region. The main pulsation period was 
found to be 363 s at a relatively large amplitude of 29 mmag. A 
second period of 382 s has been shown to be present at a much 
smaller amplitude of 3.8 mmag. This paper reports on the dis- 
covery of additional variations at a longer time scale of ~ 1 h, 
and addresses in particular the implications resulting in the case 
of an interpretation of the period as a stellar oscillation. 

2. Discovery and follow-up observations 

Photometric observations of HS 0702+6043 were obtained on 
three consecutive nights in December 1999 using the Calar 
Alto 1.2 m telescope. Details of the set-up, the observing log 
and the data reduction may be found in DR02. The light curve 
has an average cycle time of 18 s and was taken in white light. 
The analysis of the relative light curve obtained from aperture 
photometry yielded two periods in the range of 360 -390 s as 
detailed in the previous section. 

While testing a new set of time-series analysis tools (Huber 
et al., in preparation), we discovered a ~lh variation in the 
residuals of the discovery light curve of HS 0702+6043. These 
data are however of insufficient quality to conduct a false 
alarm probability analysis. In order to confirm the variation, 
we reobserved the star during two nights using the Steward 
Observatory 2.3 m Bok telescope. 

These observations, obtained in February 2004, span 
roughly 6 h on both nights. A relative light curve was extracted 
from images that have an average cycle time of 86 s and were 
taken through a F555W filter. As in DR02, the two shorter pe- 



3.1. The binary hypothesis 

From the uncontaminated optical spectrum of HS 0702+6043, 
we can easily exclude a potential F- or G-type main sequence 
companion. Because the star is relatively faint, 2MASS photom- 
etry is only accurate enough to require that any hypothetical 
companion must be later than K-type. At the supposed 1 h or- 
bit, the Roche geometry marginally allows to fit an object near 
the minimum of the mass-radius relation, i.e. an old massive 
brown dwarf. However, all systems containing an sdB plus low- 
mass companion show a substantial reflection effect. The am- 
plitude of the observed variation is quite low when compared 
with known sdB plus M-dwarf close double stars. HW Vir was 
the first such system dis covered, with pea k-to-peak amplitudes 
of ~0.26mag in V ("e.g. lKiss et allEoOOh . as well as eclipses, 
whil e the most re cent discovery is that of HS 2333+3927 
llHeberetall E004). a system with peak-to-peak amplitudes of 
0.33 mag in R. Another system, PG 1017 -086, shows a lowe r 
amplitude of 0.083 and a period of 1 .75 h JjMaxted et alll2002i) . 
however this variation is still ~22 times larger than the one we 
see in HS 0702+6043. One possible way to explain a very low- 
amplitude reflection is by assuming a low inclination angle, i.e. 
so that the system is virtually seen pole on; this is also the only 
way to prevent eclipses in very tight configurations. 

If the 1 h variation we see in HS 0702+6043 is due to 
an ellipsoidal distortion of the star by a compact companion, 
then we would expect that the actual orbital period of the 
system is twice that of the intensity variation. There are two 
known sdB plus white d warf systems sho wing an ellipsoidal ef- 
fect. KPD 0422+5421 JKoen et all) 19981) and KPD 1930+2752 
(Bille res etafl hoOO: Maxted et al., 2000). These stars show in- 
tensity variations of 0.0148 mag and 0.013 mag, respectively, 
and both are seen almost edge on (with inclinations ^80°). 
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Frequency (//Hz) 

Fig. 2. Discrete Fourier transform of the full light curve, with the window function for the two-night data set displayed in the 
inset panel in the upper right corner (frequency in pHz as in the main plot, the amplitude scaling is arbitrary). fa and fa clearly 
stand out slightly below 3000 //Hz, as does the unresolved fa conglomeration near 300 //Hz. Note the "shoulder" of that feature 
towards higher frequencies which is responsible for significant residual power after prewhitening of one low-frequency period. 



Feige 48, an EC 14026 star with parameters very similar to 
HS 0702+6043, was found to be a sdB plus whi te dwarf sy stem 
of low inclination with a period of 9 h JO'Toole et allE004l) . If 
such a system were reduced to a 2 h orbital period, ellipsoidal 
variations might become apparent. In addition to the inclina- 
tion, the size of the variation depends on how much of its Roche 
lobe the sdB fills. These relaxed constraints on the amplitude 
together with a larger possible orbit (due to the twice as long or- 
bital period) to accommodate a second component that is more 
compact make it harder to rule out the sdB plus white dwarf 
scenario. 

However, if we examine the light curve after the two higher 
pulsational frequencies have been subtracted (the thick line in 
Fig.^, it becomes quite evident that the amplitude of the varia- 
tion is changing. This strongly suggests that the variation is not 
due to binary motion. 

3.2. A rotating star spot? 

Another possible cause to be considered as the origin of the 
variation in HS 0702+6043 is a rotating star spot. Although sin- 
gle sdBs do not generally rotate very rapidl y, with only one pos - 
sible exception (PG 1605+072, but see iKuassivi et alJ Eo05). 
the following scenario can be constructed around this possibil- 
ity. Assuming that the 1 h period corresponds to the rotational 
period of the star, the variation would then be caused by an 
inhomogeneous surface temperature distribution (possibly re- 
sulting from a magnetic field). The variation has been seen at 
the same period in two observations taken four years apart, a 
time scale that is probably much longer than the typical life 
time of a star spot. However, newly appearing star spots are 
bound to result in the same rotation period determination, only 
with arbitrary phase discontinuities. Furthermore, differential 
rotation could explain the variation in the period ("multiperi- 
odicity") when a single spot drifts, or when several spots are 
present simultaneously. 

For the frequencies in the p-mode domain, substantial ro- 
tational splitting due to the fast rotation would have to be taken 
into account for all Z + modes. For a hypothetical / = 1 mode 
/o at a suitable geometric orientation, the two m — + 1 frequen- 



cies could be observable and would lie at about fa = fa + Af 
and fa = fa - Af, where Af would correspond to the rotational 
frequency. Using f and fa from Tableland fa = [fa + fa) /2, 
Af would amount to 74 fiHz, only a fraction of fa . 

The fact that no frequency splitting at the large value of 
fa (or multiples of it) is seen does not automatically rule out 
this interpretation, but, together with the unusually high rota- 
tion rate implied, makes it relatively unlikely. 

3.3. The long-period pulsation hypothesis 

The final possible cause of the observed variation is pulsation. 
iGreen etalJ §003) first detected oscillations in several sdBs 
with periods from around 30 to 80 minutes. The periodicity we 
observe lies within this range, and the amplitude is also very 
comparable to that of the IpsdBVs. As noted in the previous 
sections, the data imply multiperiodicity, a strong indicator for 
a pulsational origin of the variation. This might already be seen 
in the light curve, especially in the top panel of Fig.^ where 
the amplitude in the second half of the curve appears to reduce 
almost to zero - suggesting beating between several modes. 
The unresolved multiperiodicity becomes particularly evident 
after the highest peak in the low-frequency domain (fa) has 
been removed: After this additional prewhitening of the strong 
variation, there is significant residual power remaining, above 
the 3 cr level derived from a false alarm probability analysis. 

It is also interesting that the low frequency we measure is 
very close to, but not exactly, the difference between the two 
higher frequencies presented in Tabled This near-resonance is 
striking, and again indicates that the variation is most probably 
not due to rotation or orbital motion. If either of these were the 
case, we would expect exact resonance. All of this together sug- 
gests the very exciting possibility that HS 0702+6043 shows 
both p-modes and g-modes simultaneously. 

4. Implications for asteroseismology 

Stars that show observable multimode pulsations in both p- 
and g-modes have an enormous potential for asteroseismology 
since modes on the acoustic and gravity branches, respectively, 
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probe different regions in the depth-dependent stellar structure. 
HS 0702+6043 is probably one among very few objects to fall 
into this class of hybrids. Not only does the range of its long 
period(s) and the magnitude of the corresponding amplitude 
match what is expected for g-mode pulsators, its stellar param- 
eters also place it at a position at the edge of both of the empir- 
ical instability regions in the T e fr - logg diagram (see Fig.[3j. 
In fact, from its position in Fig.|5] which shows the two groups 
of hot, short-period pulsators, and the cooler, long-period pul- 
sators, one could argue that g-mode pulsations should almost 
be expected in HS 0702+6043 (although it should not be for- 
gotten that both types of pulsators among sdB stars co-exist 
with a larger fraction of non-pulsating sdBs in the same areas). 

A second candidate for the class of sdBV/lpsdBV hybrids 
established through the initial announcement of o ur discovery 
JSchuh et alll2005b has rece ntly been published b y fearan et all 
]2005b : Balloon 090100001 JOreiro et aUl20()4 . It is striking 
that there appears to be a gap along the EHB, so that the four 
coolest sdBVs (Feige 48, HS 2201+2610, Balloon 090100001 
and HS 0702+6043) are closely connected in parameter space 
to the IpsdB Vs rather than to the rest of the sdBVs. In fact, two 
of them are hybrid (HS 0702+6043 and Balloon 090100001), 
so what about Feige 48 and HS 2201+2610? If attempts to de- 
tect long-period variations remain unsucessful, the latter two 
stars would belong to the sdBs in the IpsdBV regions that do 
not show g-mode pulsations. From a theoretical point of view, 
it is a challenge to drive both p- and g-modes simultaneously in 
this object; future models will have to account for the existence 
ofHS 0702+6043. 

On the main sequence, where the sdBV/ IpsdBV groups 
might have an analogy in the /3 Cep / [SPB] var iables, poten- 
tial hy brid s tars include two objec ts presented bv ljerzvkiewiczl 
1 1993b and lHandler et all J2004). The c lass o f yDor/dScu 
hybrids are discussed in | Handler et alJ (2002), Dupr et et alJ 
fc005b . and lHenrv & FekelN2004 . 

5. The future 

The amplitude variability observed in HS 0702+6043 is highly 
suggestive that the 1 h variation is not due to binary or rota- 
tional motion, and results from g-mode pulsations instead. 

To resolve the suspected several low-frequency modes of 
HS 0702+6043, longer observations are clearly required. To 
this end, we have recently completed a 10-day-multisite cam- 
paign to observe HS 0702+6043 photometrically that we will 
report on in a forthcoming paper. 
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